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.Abstract 

o 

[ We first derive the mass sum rules for the heavy hybrid mesons to obtain the 

binding energy and decay constants in the leading order of heavy quark effective 
theory. The pionic couplings between the lightest 1 h hybrid (Qqg) and the lowest 
three heavy meson doublets are calculated with the light cone QCD sum rules. With 
SUf(3) flavor symmetry we calculate the widths for all the possible two-body decay 
processes with a Goldstone boson in the final state. The total width of the 1 
hybrid is estimated to be around 300 MeV. We find that the dominant decay mode 
of the 1 h heavy hybrid is 1 — > 7r + 1 + where the 1 + heavy meson belongs to the 
(1 + , 2 + ) heavy meson doublet. Its branching ratio is about 80% so this mode can be 
used for the experimental search of the lowest heavy hybrid meson. 
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^ 1. Introduction 

From quark model we know that a qq meson with orbital angular momentum 1 and to- 
tal spin s must have P = (— 1)' +1 and C = (— l) l+s . Thus a resonance with J PC = 
, + ~, 1~ + , 2 + ~, • • • must be exotic. Such a state could be a gluonic excitation such as 
hybrids, glueballs or multiquark states. The hybrid and glueball has been a missing link 
in the hadron spectrum. Recently there appears experimental evidence for a J PC = 1 h 
exotic JTJ, 0, |3|]. The emergence of evidence for hybrids indicates the presence of dynamical 
glue in QCD and will be a direct test of low energy sector of QCD. 

The light 1~ + hybrid meson mass has been found to lie between lGeV ~ 1.7GeV [0-|| 
with QCD sum rules 0. With the same method the hybrid meson containing one heavy 
quark and heavy hybrid quarkonium has been studied in ||. Some decay modes of the 
light 1 h hybrid meson are discussed in |J. Recently the masses of the hybrid quarkonium 
(QQg) states are calculated in the limit of mg — > oo [[UJ. 

The hybrid meson spectrum have been studied extensively with other theoretical ap- 
proaches including the constituent gluon models ||11|| , the flux tube models 0], the MIT 



bag model |[L3| , and the lattice gauge theory yjfl. Their decays have been studied in the 
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constituent gluon model via constituent gluon dissociation [TJj , in the different versions of 



the flux tube model with a 3 Pq pair creation mechanism [TE\. The constituent gluon model 



predicts a significant width for the r/ii channel for the light 1 + hybrid meson [15| while 



the width of this channel is very small in other models H, In contrast the flux tube 
model predicts very characteristic decay modes for the 1 h decays. The dominant modes 
are 1 h — > bxir, fan or a^. In this model the width of the pir channel is also significant. 

The QCD sum rule was used to calculate the pir, i]tt channels of the 1 h hybrid mesons 
with the three point correlation functions at the symmetric point ||. The final sum rule 
for the decay coupling constants suffers from large continuum and excited states contami- 
nation. 

Since most of the decay calculation has focused on the light hybrid mesons, it will prove 
valuable to calculate the mass and decays of the hybrid meson with a heavy quark to see 
whether the same characteristic decay modes exist as in the flux tube model. 

The combination of the heavy quark effective theory (HQET) [TJJ and the light cone 
QCD sum rules Jl8| provides a convenient framework to calculate the strong and electro- 



magnetic decays of heavy hadrons containing one heavy quark [19|, [2C], |21 . 

In this work we first derive the mass sum rules for the heavy hybrid mesons and perform 
the numerical analysis. Then we calculate the pionic couplings between the lightest 1 h 
hybrid (Qqg) and the lowest three heavy meson doublets. Invoking SUf(3) flavor symmetry 
for the coupling constants we calculate the widths for all the possible two-body strong decay 
processes with a Goldstone boson in the final state. 

2. The mass sum rules for the heavy hybrid mesons 
in HQET 

2.1 Heavy quark effective theory 

The effective Lagrangian of the HQET, up to order 1/m.Q, is 

£ eff = h v iv-Dh v + — *- K + S + 0{l/m 2 Q ) , (2.1) 
2m Q 2m Q 

where h v (x) is the velocity-dependent field related to the original heavy-quark field Q(x) 
by 

h v (x) = e m ^ x ^^Q(x) , (2.2) 

f M is the heavy hadron velocity. JC is the kinetic operator defined as 

K = h v {iD t fh v , (2.3) 

where D£ = — (v ■ D) v^, with D 11 = <9 M — ig is the gauge- covariant derivative, and 
S is the chromomagnetic operator 

5 = | C mag {m Q /n) h v o^G^K , (2.4) 



' / \ \ 3/A) 
a s (fi) 



where C mag = ^-^T > A) = H - 2n//3. 
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2.2 The derivation of the sum rules 

The interpolating current for the J PC = 1 h , ++ heavy hybrid mesons in HQET reads 

J,{x) = q(x)g s YG%(x)^-h v (x) , (2.5) 



and for the J PC = 1 + ,0 ones 



= q(x)g s Yl,G%(x)-h v (x) . (2.6) 

We consider the correlators 

% J d 4 xe* kx (0\T{J^(x), 4(0)}|0) = -{g„ v - v^)Ui(u) + v,v u U 2 (u) , (2.7) 

i J d 4 xe lkx (0\T{ Jj(x), J, 5t (0)}|0) = -{g, v - v,v u )U 3 (co) + Vfl v v U^u) (2.8) 

with io = k ■ v. The imaginary parts of rii,n2,n3,n4 receive contributions from the 
J = 1 h ,0 ++ , l + ~,0 hybrid intermediate states respectively. 

To simplify the notations we denote the J PC = 1 h , ++ , l + ~, hybrid mesons with 
up or down quark by Hi, H 2 , H 3 , if 4 and those with strange quark by H*, H%, HI re- 
spectively. We define the overlapping amplitude fi as 

(01^(0)1^) = /^, (2.9) 

(0\J,(0)\H 2 ) = f 2 v, , (2.10) 

(0|jJ(0)|fl3>=/3^ , (2.H) 

<0|^(0)|iJ 4 ) =ihv,, (2.12) 

where e 1 , is the Hi, H 3 polarization vectors. 
The dispersion relation for Hi(u) reads 

YlJu) = / ds , (2.13) 

J s — ijj — le 

where p(s) is the spectral density in the limit mq — ► oo. 
At the phenomenological side 

1 ft 

Hi(u) = -- — h excited states + continuum . (2-14) 

2 Aj — to 

In order to suppress the continuum and higher excited states contribution we make 
Borel transformation with the variable uj to ( |2.13D . We have 



1 n A, f s a 

~fie-- = / P*(s)e-Tds , (2.15) 

2 Jo 

where Sq is the continuum threshold. Starting from s we have modeled the phenomeno- 
logical spectral density with the parton-like ones including both the perturbative term and 
various condensates. 
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We need the spectral density Pi(s) at the quark level. The relevant feynman diagrams 
for the derivation of Pi(s) are depicted in FIG. 1. The first line is the heavy quark prop- 
agator. The broken solid line, broken curly line and a broken solid line with a curly line 
attached in the middle stands for the quark condensate, gluon condensate and quark gluon 
mixed condensate respectively. We consider condensates with dimension less than seven. 
The last four diagrams in FIG. 1 involve with quark gluon mixed condensates. They appear 
as a s (qg s <j ■ Gq). Compared with the quark condensate, they are typically suppressed by a 
factor where ??1q = ^ 9s ^ q ^ = 0.8GeV 2 and T ~ lGeV. Its contribution is negligible. 
At dimension six there are two condensates, the four quark condensate and triple gluon 
condensate, corresponding to the sixth and fourth diagram in FIG. 1. But in the present 
case the four quark condensate appears as a 2 s (qq) 2 . It is of high order in ot s so it can be 
omitted safely. Recall in the QCD sum rule analysis of the light hybrid meson masses, the 
four quark condensate plays a cruicial role |J. Up to now there exists no reliable way to 
estimate its value. Two approaches are commonly used to deal with this problem. One is 
to invoke the vacuum saturation hypothesis and use the factorization approximation. 
The other is to scale the value derived from the vacuum saturation hypothesis by a number. 
In a factor two is used to estimate the four quark condensate value, which introduces 
large uncertainty. In contrast, in the framework of HQET, the dominant nonperturba- 
tive corrections in the QSR analysis of heavy hybrid meson masses are due to the quark 
condensate and gluon condensate, which have been determined rather precisely. 

In the calculation we need the following formulas for the gluon condensates. 

< 9 2 8 Gap G lv >= ^r(9a^9/3u - 9au9^) < 9 2 S G 2 > , (2.16) 



< 9\i ahc G a iliV G h Ci f i G c p(J >= — < glf abc Gt jS G b Se G c n > (g^9au9fip + 9p,/39ap9au 

+gaa9w9v{3 + 9pugfj,a9^(T ~ 9 p/39 'aa 9 ' pv ~ 9^9ap9vii ~ 9av9pp9(ia ~ gpp9p,a9ua) ■ (2.17) 

(\a \b \c\ 1 
T ~J = Z (<U + </-*), (2.18) 

where d a b c , f a bc are the symmetric and anti-symmetric SU(3) color group structure con- 
stants. 

The heavy quark propogator has a simple form in coordinate space. 

_ POO 

<0\T{hJx),hJ0)}\0>= dtS(x-vt). (2.19) 

Jo 

It's convenient to calculate the Feynman diagrams directly in coordinate space. Then we 
perform Wick rotation and make Borel transformation with the variable u using the Borel 
transformation formula 

Ble^ = 5(a - i) . (2.20) 

As a last step we make a second Borel transformation to IT(T) with respect to r = 7^ to 
get the spectral density, 

#00=2^11^=1/^. (2.21) 
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Finally we get 



+C4 I ^k 3 G 3 ) + c 5T ^(^ 2 G 2 ) S 5( S ) , (2.22) 

where a q = -4:ir 2 (qq), q = u,d,s and (g 3 G 3 ) = (g 3 fabcG^G^G^) . The coefficients c* are 
collected in TABLE I for 

The contribution from the nonzero light quark mass is calculated up to the order 0(m q ). 
This term is negligible for sum rules involved with up and down quark but it is important 
when the heavy hybrid meson contains a strange quark as we shall see later. 

For the leading order binding energy of heavy hybrid meson we have 



Jq" sp(s)e rds 
fo° p(s)e~*ds 



2.3 Numerical analysis of the mass sum rules 

We use a s = = — — —72 and kept four active flavors with Anrn = 220 MeV. The 

(ll-fn f )ln( . 80/2 ) 2 ^ 
3 1 QCD 

tiny up and down quark mass is taken to be zero. The quark and gluon condensates adpots 
the standard values 

(uu) = -(0.225 GeV) 3 , 
(a s G 2 ) = 0.038 GeV 4 . (2.24) 

The value of the triple gluon condensate is not well known. In fact several values exist in 
literature. We use |7[ 

(g^G 3 ) = (1.2 GeV 2 )(a s G 2 ) = 0.045GeV 6 , (2.25) 



which is smaller than the value 0.06 — 0.1 GeV 6 in |22 and a even larger value 0.4 GeV 6 
from the "instanton liquid" approach to the QCD vacuum Q. Later we will enlarge ( [2.25D 
by a factor of ten to see the uncertainty from this source. 

There are two commonly used methods to extract the masses, the derivative method 
and the fitting method. With the derivative method we arrive at ( j2.23|) . The fitting method 
involves with fitting the left hand side (L.H.S.) and right hand side (R.H.S.) of Eq. ( p,15|) 



with the most suitable parameters Aj, fi, s° directly in the working region of the Borel 
parameter. In the numerical analysis, we invoke both methods to crosscheck our results. 
We find both methods yield nearly the same results. 

We require that (1) the absolute value of each condensate contribution be less than 30% 
of the leading perturbative term with continuum subtracted and (2) the sum of the power 
corrections be less than one third of the whole sum rule. This requirement leads to the 
lower limit of the continuum threshold. Typically the triple gluon condensate is less than 
0.5% of the leading term with the value in ( j2.25| ). The series of the operator expansion 
converge fast. 

We present the numerical values of A, /, sq in TABLE II. It is understood that there is 
an error of 0.1 GeV for A, s - With these values the left hand side and right hand side of 
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2.15| ) agree within one percent in a large interval of T. Correspondingly, the ratio between 



the quark, gluon, triple gluon condensate and the perturbative term is collected in TABLE 
III. 

The dependence of the binding energy on the Borel parameter T with different contin- 
uum threshold So is shown in FIG. 2 and 3 for the two exotic hybrid mesons. We have also 
plotted the fitting lines and the curves of the right hand side of Eq. (|2.15 ) for Hi, H4 in 
FIG. 4 and 5. 

Varying Aqqq from 220MeV to 300MeV the final result changes within 5%. in our 
numerical analysis. 

If we use {g 3 G 3 ) = 0.45GeV 6 , the values of A, so is shifted upwards by 2% for Hi. For 
the other three channels the best fitting parameters changes little. In other words, the 
uncertainty due to the triple gluon condensate is small and included in the errors given 
already. This is in strong contrast with the case for the light hybrid mesons, where the 
hybrid masses depend cruicially on the value of the triple gluon condensate. Varying its 
value from 0.045GeV 6 to 0.4GeV 6 , the 1 h hybrid meson mass increases from 1.0 GeV to 
1.5 GeV ||§. 

For the numerical analysis of the strange heavy hybrid meson, we use (ss) = 0.8 (mm), 
m s = 150 MeV. The results are collected in TABLE IV and V. As can be seen from TABLE 
V the strange quark mass correction is very important. The binding energy of the light 
component of Hf is roughly 130 MeV larger than that for Hp The continuum threshold 
increases by about 150 MeV. The decay constants increase typically by 25% due to the 
strange quark mass correction. 

It is interesting to note that the central value of Ah is much greater than that for the 
(0~, 1~) doublet of Qq meson, A_ 1 = (0.5 ± 0.1)GeV. And the continuum starts at rather 
large s due to the presence of the dynamical gluon in Hq. 

The masses of heavy hybrid mesons shall be around (m& + Ah). If we can derive Ah 
reliably, we have a good estimate of Hq mass. Especially for the bottom quark system, 
the 1/ttlq correction is not large. This point has been the motivation of our considering 
heavy hybrid mesons with one heavy quark in the framework of HQET. 

The uncertainty due to the dimesion six condensates renders the reliable extraction 
of the light hybrid meson masses rather difficult. In this section we have calculated the 
binding energy of the heavy hybrid meson masses in the leading order of HQET with 
QCD sum rules. Within the present approach (1) the heavy quark mass is disentangled; 
(2) the four quark condensate is of higher order in a s , hence can be neglected safely; (3) 
the dominant power corrections are from the quark condensate and the gluon condensate, 
which is well known. The triple gluon condensate shall at most affect the binding energy 
by 3% varying its value from 0.045GeV 6 to 0.4GeV 6 . 



3. Light cone QCD sum rules for the pionic couplings 

Note Hi is the lightest exotic heavy hybrid meson while H± lies about 1.4 GeV higher 
than Hi. Experimental discovery of the ++ and l + ~ heavy hybrid mesons will be difficult 
since they have the same non-exotic quantum numbers as the radial excitations of ordinary 
heavy mesons. In this section we discuss the decay modes and widths of the loweset heavy 
hybrid meson Hi. 

Denote the doublet (l + ,2 + ) with j e = 3/2 by (Bi,B*), the doublet (0 + , 1+) with 
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j t = 1/2 by (B' , B[) and the doublet (0~, 1") by (B, B*). 
We need the following interpolating currents [B3[ 



= ^ f WH) (c? - ~7?^) ? , (3-1) 
4+| 2 = (tW + iT^T - \g? a% a) ? , (3.2) 

jj" 1 = fc B 7t t 9 > 4" I = K^q , (3-3) 



where 7f = 7 M — g 1 ^ = — ^"iv We use "t" to denote the tansverse index. 

We take the p-wave decay process H\ — > Bit as an example to illustrate the method. 

• H l -> Bit 

The decay amplitude is 

M{H X - Bn) = eUlgi , (3.5) 

with g* = - (q ■ 

For deriving the sum rules for the coupling constant gi we consider the correlator 

i J d 4 x e**(v{q)\T (Jo,-,i^)^(0)) |0> = G^w, u/)g< • (3.6) 

The function G\(uj,u}') in ( |3.6| ) has the following double dispersion relation 

f-ifH9i c d , 

+ 7 + ~, ~, . (3-7) 



4(A_i - a>) (Ajj- — a/) A_i - A H - u' 
where Apj £ = mp J( — mg and fpj e are constants defined as: 

(0| ^(0)1/, P',j'i) = fpiAfSpp^rj^ . (3.8) 

Keeping the three particle component of the pion wave function, the expression for 
G\{uj,uj') with tensor structure reads 

Gi(lo,u') = J d 4 xe lkx S(x - vt)Tr{ l5 ^Y < ^(g)k(0)^(0)g(x)|0 >} , (3.9) 
where we have not included the two particle component of the pion wave function 

2 

since they are of higher order in a s and suppressed by a large factor , which arises 
from the additional loop integration with the gluon attached to one of the quark line. 

The light cone three particle pion wave functions are defined as [ |19|] : 

< Tt(q)\d(x)a a pj 5 g s G IMl/ (ux)u(0)\0 >= 
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< n(q)\d(x)-f^ 5 g s G a(3 (vx)u(0)\0 >= 

+U-^-(q a x - q/sx a ) [ Va in (a i )e i ^ +V ^ (3.11) 
a • x J 



and 



< n(q)\d(x)-f^g s G af} (vx)u(0)\0 >= 

ifMgc* - ^) - qJgn - / Va ti p ± ( ai )e^ + ^ 

LV q ■ x ' v q • x n J 



+iU-^(q a x (3 - q p x a ) J Va i ip^a i )e v * B ^ +va *> . (3.12) 

The operator G a p is the dual of G a p\ G a p = 2 ^ a pspG 5p ] T>oti is defined as Da^ = 
daida 2 da 3 S(l — a\ — a 2 — 03). Due to the choice of the gauge x^A^x) = 0, the 
path-ordered gauge factor Pexp (ig s fo dux^A^ux)) has been omitted. 

The function y? 37r is of twist three, while all the wave functions appearing in eqs. (|3.11|) , 
Q3.12| ) are of twist four. The wave functions ip(xi, ii) (li is the renormalization point) 



describe the distribution in longitudinal momenta inside the pion, the parameters Xj 
{J2i x i = 1) representing the fractions of the longitudinal momentum carried by the 
quark, the antiquark and gluon. 

The wave function normalizations immediately follow from the definitions / Da^^iai) 
1, fVai<f±(ai) = J Vatican) = 0, JVai^ ± (ai) = - / ' Va^^ai) = 5 2 /3, with the 
parameter 5 defined by the matrix element: < Ti{q)\dg s G a ^ a u\Q >= i5 2 f^q^. 

Expressing (|3.9|) with the pion wave functions we arrive at: 



Gti^Lo') = -J dtj P aj e^ + ^)H/ 3 ^(a^(g^) + /.[^(a l )-^p]} + 

(3.13) 

For large Euclidean values of u and uo 1 this integral is dominated by the region of 
small t, therefore it can be approximated by the first a few terms. 

After Wick rotations and making double Borel transformation with the variables u 
and u' the single-pole terms in ( |3.7| ) are eliminated. We arrive at: 



IgiL^fHe-^^ = fA{uo)T 2 - A-[*jl(«o) " ^^]T , (3.14) 

where /# — fi, Uq — T J ^ T , T = , 7i, T 2 are the Borel parameters. Note the sum 
rule is asymmetric with the Borel parameter T\ and T 2 . The continuum subtraction 
is complicated in the present case. We shall discuss this point in section |5]. The right 
hand side of the sum rule is the result after integration of the double spectral density 
with respect to s±, s 2 in the interval (0, 00) so it includes the continuum contribution. 

The new wave functions introduced in (|3.14j ) are defined as 

$i(ax) = (p i (a 1 ,a 2 ,l - ax - a 2 )da 2 , (3.15) 

Jo 
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with i = 3tt, _L, || etc and 



$3^( M o) = — ^— — |a 1= « • I 3 - 16 ) 



We have used the Borel transformation formula: B^e auJ = 5(a — ^). Integration by 
parts are employed to absorb the factor (q-v), which leads to the derivative in ( |3.14| ). 
In this way we arrive at the simple form after double Borel transformation. 

M{H X - B*n) = ie IMX(J pete a q^g 2 , (3.17) 
where is the polarization vector of B*. 
Similarly we consider the correlator 

i f d A x e ik - x (-K(q)\T (j^i (x)J^(O)) |0> = i W gVG 2 ( W , W ') . (3.18) 



G 2 (u,u') = k jHdi J Va i e« kv * IW * 1 '> t [<p ± (a i ) - <M«i) + <p\\(ai)] +■■■ . (3.19) 

1 A - 1 A f 

-^/^//fe"^^' = -y[<MO - $±M + $||K)]T . (3.20) 

• # x -> 

This process is forbidden due to the parity and angular momentum conservation. 

There exist two independent coupling constants, corresponding to S-wave and D-wave 
decay. The decay amplitudes are: 

M(Hx - = if (£{^03 + - gfrJJM . (3.21) 

where t]^ is the polarization vector of B*. 
We consider the correlator 

i J d'x e^(n(q)\T(jl +l2 (x)J^(0)) |0) = gfcGfaoS) + (£ g J - rfflJJG^.o/) . 

(3.22) 

G? s (w,w / ) = /J^/^^^ . 

(3.23) 

G 4 (u;,u/) = [°° dt [ Va^^ {h^a^q ■ v) + \J^-[ l p L {a i ) 
Jo J 2q ■ v 

+<p\\(ai) + ^±(«i) + ^n(oi)]} + ■ • • • (3.24) 
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i S 3/ + ,j/»e-<^» = - hl ^^ {m ^- m( u o) + ^r ■ (3.25) 



where 



^4/ + ,i/^e- ( ^ + ^ ) = h^{u )T + ^(u ) , (3.26) 
$(«) = J da[^ ± (a) + ^||(a) + ^ ± (a) + ^|(a)] . (3.27) 



#1 -> i?i7T 

There also exist two independent coupling constants, corresponding to S-wave and 
D-wave decay. The decay amplitudes are: 



M{H, -> B l7r ) = I7°cf {^flfe + - qWM ■ (3-28) 



where e M is the polarization vector of B\. 
We consider the correlator 



i J d'x e ik -*{n(q)\T(j° + s(x)J^(0)) |0> = g^G.M + - <&k)G 6 (a;,u/) . 

(3.29) 



G 5 (uj,uj') = -^yj^dt J Va^+^a^q-vfih^a^q-v) 

+fA<P±(ui)- n ^}} + -- - ■ (3-30) 



G 6 (cj,u/) = ^ j™dt J ^(^"^{^feKKg^l + A^ila,) 

-2^||(ai) - <p±(aii) + <p\\(aii)]} H . (3.31) 



^5/+,|/ffe T i T 2 = — {/W ^3 L=« T 



^6/ + ,3^e ^ * ="^-{4/3. ^ U= U0 T 

-A«o[5$±(«o) - 2$||K) - $l(u ) + $|| W]T} . (3.33) 
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Hi -> B*% 

There exists only one independent coupling constant, corresponding to D-wave decay. 
The decay amplitude is: 

M(#! -> £? 2 *tt) = ^eJfi^Xg^ + («i <-> a 2 )]<? 7 • (3-34) 

where ?7 aia2 is the polarization tensor of 
We consider the correlator 

i J d A x e lk - x (n(q)\T (j 2 a +| (x) J^(0)) |0) = [ie* 1 ^^ + («i <- a 2 )]G 7 (w,£</) . 

(3.35) 



G 7 (u,uj') = ^ r dt /Pa J e i ^+^ 1 Vb±(a i )-^(a J )+^||(a. J )] + --- . (3.36) 



4 7o 





(uo)]T . (3.37) 


4. Determination of the parameters 




4.1 The values of A, / 




We need the mass parameters A's and the coupling constants /'s 
interpolating currents as input. The results are |[24|. p5j| 


i of the corresponding 


A_ i = (0.5 ± 0.10)GeV , 


(4.1) 


f_i = (0.35 ± 0.04)GeVi , 


(4.2) 


A +i i = (0.85±0.10)GeV, 


(4.3) 


f + x = (0.36 ± 0.04)GeVi , 


(4.4) 


A +> 3 = (0.95±0.10)GeV, 


(4.5) 


/ +) 3 = (0.28 ± 0.03)GeVt . 


(4.6) 


The value of the continuum threshold for the (0~, (0 + , 1 + 
(1.1 ± 0.1), (1.2 ± 0.1), (1.3 ± 0.1) GeV respectively. 


), (1+2+) doublets is 


4.2 Expressions of PWFs 




The detailed expressions of the pion wave functions relevant in our 


calculation are: 


y?37r(ai) = 360«ia 2 a3[l + ^(7a 3 - 3)] , 


(4.7) 


<p ± {oi) = 305 2 ( ai - a 2 )a 3 2 [i + 2e(l - 2or 3 )] , 


(4.8) 
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,;) = 1205 2 e(a 1 — a 2 )a;ia; 2 a3 , 
305 2 (l-a 3 )a 2 [i + 2e(l-2a 3 ) 



e(l - 2a 



3JJ , 



(4.9) 
(4.10) 

(4.11) 



where the coefficients ^3^, u 10 [PE], 5 2 [27| and e [|28| have been determined from QCD sum 
rules. At the scale fi = 1.0 GeV, cj 1>0 = -2.88, / 37r = 0.0035GeV 2 , 5 2 = 0.2GeV 2 , e = 0.5. 

With (|4.7|) -( CTD we can calculate the explicit expressions of the other wave functions 
defined in this work. Note only the asymptotic form of these wave functions are exactly 
known. The pieces involved with e, arise from nonperturbative corrections. They are 
estimated from the moments of the pion wave functions with QCD sum rules. For our 
purpose it's enough to keep the asymptotic form. 



$3*(u) = 30w(1-w) 4 , 


(4.12) 


[$3.(m)] / = 30(1-5 M )(1- M ) 3 , 


(4.13) 


(u)Y = 60ufl — 3u)(l — u) 3 


(A 14) 


[^3n{U)\ — — — DU){1 — U) , 


(A 1 ^ 


[u$ 3n (u)]'" = 720(u - l)(5u 2 - 5u + 1) , 


(4.16) 


$±H = ^ 2 (5«-l)(l-«) 3 , 


(4.17) 


u$ ± (u)}" = -^<5 2 (25w 3 - 48m 2 + 27m - 4) , 


(4.18) 


$n(n) = [$||(u)]' = [u$||(it)]" = , 


(4.19) 


$±(u) = ^6 2 (3u + l)(l-u) 3 , 



(4.20) 


[$ ± (w)]' = -10<5 2 m(1 -uf , 


(4.21) 


20 

$,|H = - y 5 2 n(l-n) 3 , 


(4.22) 


2D 

[®\\{u)]> = -^8 2 {l-4u){l-u) 2 , 


(4.23) 


$(w) = , 


(4.24) 



where '," , ' denotes the first, second and third derivative with respect to u. 
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5. Subtraction of the continuum 

We have the double dispersion relation for Gi(u,u/), 

P(si,s 2 ) f°° , Pi(si) 



jo jo [Si — u> — ze){S2 — u> — ie) Jo 



s\ — u — le) 



JO (S 2 — Co> — 26J 

where the ellipse denotes the subtraction terms. 

Making double Borel transformation to the variables u,u' we get 

/*oo roo 

n(r 1 ,r 2 )=/ ds x ds 2 e- SlT1 e- S2T2 p( Sl ,s 2 ) . (5.2) 

JO JO 

The single pole and subtraction terms have been eliminated in (|5.2[) . Making a second 
double Borel transformation to r 1 ,r 2 , we get the double spectral density 

P (s 1 ,s 2 )=B%B%IL(t 1 ,t 2 ). (5.3) 

After expressing Gi(u,u') with the PWFs and finishing Wick rotations and double 
Borel transformation, we can get the following general formula: 

/•oo rl 

n(n, r 2 ) = dt S[(l - u)t - Tl )S(ut - r 2 )t- n ip(u) 

JO JO 

T 2 x 



= ll>(- 

(T 1 +T 2 ) n+lY T-l + T 2 J 

oo ^fc 

= + (5 - 4) 

if we assume 

oo 

iP(u) = Y j a k (l-u) k . (5.5) 

fc=0 

In order to faciliate the numerical analysis we collect the coefficients a& in TABLE VI after 
expanding the PWFs into the polynomials of (1 — u). 

OO k-\-Tl C\ 

Introducing new variables s + = S1 1 2 S2 , s + = Sl + S2 , ^_ = ^ - ^, we have 

r°° s + r s + 1^- 

n(Ti,T 2 )=2/ ds+e~— ds^e T - p(s+,s_) , (5.7) 

JO J-s+ 

where 

^ -E^w^ { -^ m2s - ] ■ (5 - 8) 

The general quark-hadron duality holds only after the integration with respect to the 
variable s_ in ( |5.7| ). Finally we get 

n^,T0 = E^j^i^^i , (5.9) 
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where E c is the continuum threshold. From its definition we know that E c = s °t, s ° where 
So, sq is the continuum threshold for the mass sum rules of the hybrid and heavy meson 
respectively. 

Letting E c — > oo we recover ( |5.4| ). In the case of symmetric sum rules, only the term 
with q = k survives since T\ = T 2 = 2T and J- = 0. It is straightforward to get 

II(T)=^) J q Ec s n + e-^ds + . (5.10) 

Now the subtraction of the continuum contribution does not depend on the wave functions. 

In the present case the sum rules are asymmetric with Ti and T 2 . It's reasonable to let 
T\ = 2(3 Aq, T 2 = 2(3 Ah, where Aq is the binding energy of the heavy meson in the leading 
order of HQET and ft is the dimensionless scale parameter. Then we have u 

rp 2AqAh n 

- A Q +A H P- 



A 



Q 



A H +A Q ' 
We rewrite ( |5.9| ) as 

n( r) .^gle^t-rwf )} , (5.11) 

n k 

where f n (x) = 1 — e~ x p- is the factor used to subtract the continuum. Note ==- = 

fc=0 

xfrxf ~ i h \ for the ( 1+ ' 2+ )' ( 0+ ' (°~' 1_ ) doublet respectively. So only the first 
few terms in the bracket in ( |5.11|) is important. Replacing f n+q (^) by f n (^), we have 

n(r)=r . +1/ „ ( l )s - { g_J_ ( ll )9} 

^'/.(fWti,). (5.12) 

The above approximation is useful when ^ — ► 0, i.e., the decay heavy meson mass is very 
close to the parent hybrid meson mass. We shall use the exact expression (|5.11 ) in the 
numerical analysis below. 



6. Numerical results and discussions 

We now turn to the numerical evaluation of the sum rules for the coupling constants after 
the continuum is subtracted carefully. 

The variation of the coupling constants with the Borel parameter T and E c is pre- 
sented in FIG. 6-12. The curves correspond to E c = 1.4, 1.5, 1.6GeV respectively. Stability 
develops for these sum rules starting from 0.8 GeV. Numerically we have 

9xf-,\fs = -(0-12 ± 0.02)GeV 4 , (6.1) 

92Lif H = (0.078 ± 0.016)GeV 4 , (6.2) 

^3/ + ,|/// = (0.11±0.04)GeV 5 , (6.3) 

94f + if H = (0.046 ± 0.01)GeV 3 , (6.4) 
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<? 5 / +) §/h = (0-33 ± 0.08)GeV b , (6.5) 
9ef + ,sf H = (0.026 ± 0.006)GeV 4 , (6.6) 
97f+*fH = "(0.018 ± 0.004)GeV 4 , (6.7) 

where the errors refers to the variations with T and E c in this region. And the central value 
corresponds to T = 1.1 GeV, (3=1, and E c = 1.5GeV. Note the sum rules for gi,g4,ge 
changes less than 10% when either T or E c varies in the working region. Variation with T 
is about 20% for those sum rules for g 2 , g-j. Only the sum rule for g 5 has a rather strong 
dependence on the continuum threshold E c . It changes about 20% with E c = (1.5 ± 0.1) 
GeV. 

With the values of fa given in the previous sections we get 



gi = 


-(1.4±0.3)GeV- 1 , 


(6.8) 


#2 = 


(0.85±0.2)GeV" 1 , 


(6.9) 


gs = 


(1.3 ±0.3) , 


(6.10) 


g^ = 


(0.5±0.2)GeV~ 2 , 


(6.11) 


g5 = 


(5.0 ± 1.0) , 


(6.12) 


ge = 


(0.4±0.1)GeV~ 2 , 


(6.13) 


gr = 


-(0.3±0.05)GeV~ 2 . 


(6.14) 



With these coupling constants we can calculate the decay widths of heavy hybrid 
mesons. The leading order binding energy of Hi meson is A# = 1.6 GeV, which is not 
small compared with the bottom quark mass rrib = 4.7GeV. We shall take into account of 
the corrections due to the finite mj partly. In full QCD, we introduce the decay constant 
of the hybrid meson Hi as: 

(01^(0)1^) =m H F H el, (6.15) 



where F H = -jM= as wij, ^ oo. The amplitude of the decay process Hi — > Bn in full QCD 
is 

M(Hi -> Bn) = e^q^m H m B hi . (6.16) 



In the limit nib — > oo the decay constant hi in full QCD is equal to the decay constant g\ 
defined in HQET. The decay width formulas in the leading order of HQET are not useful 
in the present case. We include the finite correction and use the following decay width 
formulas. 

16-7T m H m B m H 

T(Hi - B*n) = U^^^^fglW , 
8n m H rriB*m H 

T(Hi -> B[n) = ^(^)(^^) 2 (3g 2 3 + Ag 3 g 4 \q\ 2 + 2gM A )\q\ , 
167r m H tub'Ttih 

F(Hi - Bin) = -U^)(^M 2 (3s 5 2 + Ag 5 g 6 \q\ 2 + 2g*\tfM , 
167r m H m Bl m H 

r(^i - Bin) = U^){^^)VM , (6-17) 
An tjih rn,B*mH 
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where \q\ = \J{m? H — (ijiq + m 1T ) 2 ){m 2 H — (tuq — m 7r ) 2 )/2m^, tuh, ttiq is the hybrid and 
heavy meson mass. Note summation over charged and neutral pion final states has been 
performed. In the calculation we use mg = 5.28 GeV, m B , = 5.33 GeV, m B ^ = m B +0.37 = 
5.7 GeV, rriB 1 = tub* = 5.8 GeV, m H = nib + A# = 6.3 GeV with m b = 4.7 GeV. For the 
heavy meson with strangeness, we use m Ba = m B + 0.09 = 5.37 GeV etc. 

We employ the SUf(3) flavor symmetry to relate the coupling constants for the decay 
processes Hi — > B S K, Hi — > Br] etc. For example, gu^B s K = QiiQh^b^ = We collect 
the numerical results of the decay widths for the different channels in TABLE VI. 

Summing all the two-body strong decay channels with one Goldstone boson in the final 
state, we estimate the total decay width of the lowest lying hybrid meson to be around 300 
MeV. The dominant decay mode is Hi — > Bin. Its branching ratio is about 80%. Now we 
know Bi is a narrow resonance with a width of ~ 20 MeV. So this mode can be used to 
detect the possible existence of Hi experimentally. The same characteristic decay modes 
for the light hybrid mesons have been predicted in the flux tube model [T(| . 

In summary we have calculated the binding energy of the heavy hybrid mesons in the 
leading order of HQET. With the help of the light cone QCD sum rules we have extracted 
the pionic couplings between the hybrid Hi and lowest three heavy meson doublets. Invok- 
ing the flavor SUf(3) symmetry we have estimated the strong two-body decay widths of 
Hi with one Goldstone boson in the final state. Our calculation yields very characteristic 
decay modes of Hi and confirms the earlier predictions on the particular decay modes for 
the light hybrid mesons from the flux tube model. The mixing between rjs and r]i and the 
possible contribution due to the QCD anomaly will be topics of fututre work. 
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China. 
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TABLE I. The coefficients in the expressions for the spectral density p(e). 





1 — h 
1 


u 


1 


A 

u 


Cl 


1 


3 


-1 


-3 


c 2 


1 


3 


-1 


-3 


c 3 


1 


-3 


1 


-3 


c 4 


1 


-3 


1 


-3 


c 5 


1 


-1 


-1 


1 



TABLE II. The values of A, /, s for the nonstrange heavy hybrid meson (Qqg). A, s 

7 

is unit of GeV, F is in unit of GeVa. 





1-+ 


0++ 


1+- 


o— 


A 


1.6 


2.17 


1.66 


3.0 


/ 


0.23 


0.425 


0.162 


0.98 


so 


2.0 


2.5 


1.9 


3.4 



TABLE III. The ratio between various condensates and perturbative term after the 
continuum subtraction for the nonstrange heavy hybrid meson. 





1-+ 


0++ 


1+- 


o— 


< qq > 


0.22 


0.34 


-0.26 


-0.14 


< g 2 s G* > 


0.16 


-0.23 


0.19 


-0.08 


< 9 3 S G 3 > 


0.004 


-0.004 


0.005 


-0.0006 


<qq>< g 2 s G 2 > 


0.02 


-0.004 


-0.026 


0.0006 
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TABLE IV. The values of A, /, s for the strange heavy hybrid meson (Qqg). 





1-+ 


0++ 


1+- 


o— 


A 


1.73 


2.28 


1.76 


3.15 


/ 


0.28 


0.538 


0.2 


1.04 


s 


2.13 


.63 


2.03 


3.53 



TABLE V. The ratio between various condensates (including strange quark mass cor- 
rection) and perturbative term for the strange heavy hybrid meson. 





1-+ 


0++ 


1+- 


o- 


< ss > 


0.15 


0.23 


-0.17 


-0.1 


< 9 2 S G 2 > 


0.13 


-0.19 


0.15 


-0.07 


< 9 3 S G 3 > 


0.003 


-0.003 


0.004 


-0.0005 


< ss >< g 2 s G 2 > 


0.01 


-0.003 


-0.014 


0.0004 


m s 


0.08 


0.20 


-0.09 


-0.15 



TABLE VI. The coefficients a k when the PWFs are expnaded into polynomials of (1—u). 
They are in unit of 5 2 for the twist four PWFs. 







a 2 


a 3 


04 


a 5 








-120 


150 










10 

3 


25 
6 




$j_(u) -$,_(«)+ 4>||(u) 






_20 


5 








360 


-600 










-10 


10 






$37r(«) ^ 








30 


-30 


[U$te(u)]'" 


-720 


3600 


-3600 






{u[$ ± (u) - SM]y" 


20 


-90 


250 
3 










-120 


300 


-180 


u[5$±(u) - 2$||(m) - + $|| (u)] 






20 

3 


55 
3 


35 
3 


U[$±(u) - + $||(u)] 






20 

3 


35 
3 


-5 



TABLE VII. The decay widths of different two-body decay channels for the 1~ + hybrid 
meson Hi, where combinations of Goldstone bosons and heavy mesons yield different final 
states. The unit is MeV. The minus sign means either such a decay mode is not allowed 
by the phase space or the decay width is negligible. 





B 


B* 


B'o 


B[ 


Bi 


B* 


7T 


11 


9 





25 


230 


0.05 


V 


1 


1 





1.6 






K 


5 


5 





10 
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gure Captions 



FIG. 1. The relevant feynman diagrams for the derivation of the QCD sum 
rule ( 2.15| ). The broken solid line, broken curly line and a broken solid line 
with a curly line attached in the middle stands for the quark condensate, gluon 
condensate and quark gluon mixed condensate respectively. 

FIG. 2 The variations of A with T and so for H\. From top to bottom the 
curves correspond to sq = 2.1,2.0, 1.9 GeV. T is in unit of GeV. 

FIG. 3 The variations of A with T and so = 3.5, 3.4, 3.3 GeV for H4. 



FIG. 4 The variation of the right and left hand side of Eq. ( 2.15 ) with T 
is plotted as solid and dotted curves respectively for Hi with the values of 
A, /, s in TABLE II. 

FIG. 5 The variation of the right and left hand side of Eq. ( 2. 15| ) with T for 
H 4 . 

FIG. 6 The dependence of g±f_ i/h on the Borel parameter T and the contin- 



1 -(— ) 

threshold E c after the factor T i T 2 is moved to the right hand 

side of Eq. ( g3p . 



uum 



FIG. 7 The variation of gif_ iIh with T and E c . 

> 2 

FIG. 8 The variation of 03/ , i/h with T and E c . 

"■"'a 

FIG. 9 The variation of 04/ , ifjj with T and E c . 

"""'a 

FIG. 10 The variation of 05/, afu with T and E c . 

i "' 2 

FIG. 11 The variation of gef, ifu with T and E c . 

' ' 2 

FIG. 12 The variation of 57/, a/jj with T and £? c . 
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